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Abstract 

We investigate the dust and gas distribution in the disk around HD 142527 based on ALMA 
observations of dust continuum, ^^CO J = 3 - 2, and C^®0 J = 3 - 2 emission. The disk 
shows strong azimuthal asymmetry in the dust continuum emission, while gas emission is more 
symmetric. In this paper, we investigate how gas and dust are distributed in the dust-bright 
northern part of the disk and in the dust-faint southern part. We construct two axisymmetric 
disk models. One reproduces the radial profiles of the continuum and the velocity moments 0 
and 1 of CO lines in the north and the other reproduces those in the south. We have found 
that the dust is concentrated in a narrow ring having ~ 50 AU width (in FWHM; Wd = 30 AU in 
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our parameter definition) iocated at ~ 170 - 200 AU from the centrai star. The dust particies 
are strongiy concentrated in the north. We have found that the dust surface density contrast 
between the north and south amounts to -- 70. Compared to the dust, the gas distribution is 
more extended in the radiai direction. We find that the gas component extends at ieast from 
--100 AU to ^ 250 AU from the centrai star, and there shouid aiso be tenuous gas remaining 
inside and outside of these radii. The azimuthai asymmetry of gas distribution is much smaiier 
than dust. The gas surface density differs oniy by a factor of -- 3 -10 between the north 
and south. Hence, gas-to-dust ratio strongiy depends on the iocation of the disk; ~ 30 at the 
iocation of the peak of dust distribution in the south and ~ 3 at the iocation of the peak of dust 
distribution in the north. Despite iarge uncertainties, the overaii gas-to-dust ratio is inferred to 
be 10 - 30, indicating that the gas depietion may have aiready been under way. 

Key words: stars: individual (HD 142527) — stars: pre-main-sequence — planetary systems; protoplan¬ 
etary disks — submillimeter: planetary systems — radiative transfer 


1 Introduction 

Transitional disks are circumstellar disks having an inner hole of dust emission, and are considered to be in the evolutionally phase 
from gas-rich protoplanetary disks to gas-poor debris ones (Strom et al. 1999; Calvet et al. 1999; Andrews et al. 2011). This class 
of disks have attracted much attention as valuable samples to study the disk evolution and planet formation processes. Among a 
number of transitional disk objects, HD 142527 is a subject of intense study. It is a Herbig Fe star (Waelkins et al. 1996) harboring a 
disk with a wide dust cavity with the radius of > 100 AU (Fukagawa et al. 2006; Fujiwara et al. 2006; Verhoeff et al. 2011; Rameau 
et al. 2012; Casassus et al. 2012). The stellar mass is ~ 2.2 Mq and the age is 5 Myr if we adopt d = 140 pc considering the 
association to Sco OB2 (Fukagawa et al. 2006; Verhoeff et al. 2011; Mendigutla et al. 2014). Recent observations have revealed the 
possible existence of low-mass companion (~ 0.1 — 0.4 Mq) at 13 AU from the central star (Biller et al. 2012; Rodigas et al. 2014). 

The disk around HD 142527 shows several interesting features. The near infrared scattered light image of the disk shows extended 
emission out to ^ 300 AU and large-scale spiral features are observed (Fukagawa et al. 2006; Casassus et al. 2012; Canovas et al. 
2013; Rodigas et al. 2014), hinting that some dynamical activity is taken place in the disk (Casassus et al. 2012). ALMA observations 
of dust continuum emission show significant azimuthal asymmetry with a bright horseshoe-like emission in the northern part of the 
disk (Casassus et al. 2013; Fukagawa et al. 2013; Perez et al. 2015). It is indicated that large grains are concentrated in this northern 
region (Casassus et al. 2015). In the vicinity of the central star, a stream-like features in HCO+ J = 4 — 3 emission (Casassus et al. 
2013) and a point source of dust emission (Fukagawa et al. 2013) are observed. HCN J = 4 — 3 and CS J = 7 — 6 emissions are 
also spatially resolved with ALMA (van der Plas et al. 2014). 

Fukagawa et al. (2013) presented the results of ALMA band 7 observations of dust continuum and gas emission in ^®CO J = 3 — 2 
and C^®0 J = 3 — 2. On the basis of the very bright dust continuum emission (> 20 K) in the northern part, they discussed two 
possibilities for the disk gas distribution. One is that the gas-to-dust mass ratio is less than 100 at least in the northern part and 
significant dust concentration occurs. The other is that the gas-to-dust mass ratio is standard value of 100 so the disk gas mass can 
be high enough for the onset of gravitational instability. In either case, the disk is likely to be in the process of ongoing planet 
formation. 

It is important to pin down the amount of gas and dust by means of detailed modeling. In this paper, we derive the gas and dust 
distribution based on ALMA Cycle 0 Band 7 observations of HD 142527 using more detailed modeling of gas and dust emission 
than Fukagawa et al. (2013). We derive the dust distribution from the continuum emission and the gas distribution from the CO 
emission in order to obtain the gas-to-dust ratio in the northern and the southern part of the disk. 

The paper is constructed as follows. In section 2, we summarize the results of observations. In section 3, we describe the method 
of modeling. In sections 4 and 5, we present the modeling results of dust and gas distribution, respectively. In section 6, we discuss 
implications of our best-fit model. Section 7 is for summary. 
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Table 1. Gaussian parameters (equation (1)) that fit the radiai profiies of dust continuum emission in the brightest and the faintest 
directions. The parameters ro and w are given in the unit of AU assuming that the distance to HD 142527 is 140 pc. 


PA 

Ip [Jy/asec^] 

'^Ojobs [AU] 

[AU] 

11-31° (North profile) 

1.2 

152 

51 

211 — 231° (South profile) 

0.050 

174 

51 


2 Observation 

2.1 Observation and Data Reduction 

ALMA Band 7 Cycle 0 observations of HD 142527 (RA=J2000 15h56m41.9, DEC=J2000 -42dl9m23.3) were carried out in the 
extended array configuration with a maximum baseline of about 480 m. The observations consisted of six scheduling blocks during 
the period from June to August 2012. The correlator was configured to store dual polarizations in four separate spectral windows 
with 469 MHz of bandwidth and 3840 channels each, and their central frequencies are 330.588, 329.331, 342.883 and 342.400 GHz, 
respectively, to target the molecular lines of J = 3 — 2 and C^®0 J = 3 — 2. The resultant channel spacing for the lines was 
122 kHz, corresponding to 0.12 km s“^ in velocity at these frequencies, but the effective spectral resolution was lower by a factor 
of around 2 (~0.2 km s“^) because of Hanning smoothing. The continuum data from all the spectral windows were aggregated and 
treated as a single data set of 336 GHz in central frequency and 1.8 GHz in bandwidth. The on-source integration after flagging 
aberrant data was 3.0 h. 

Calibration and reduction of the data were made with the Common Astronomy Software Applications (CASA) version 3.4, in 
almost the same way as that in Fukagawa et al. (2013). Self-calibration was performed for the continuum to improve the sensitivity 
and image fidelity, and the final gain solution was also applied to ^®CO and C^®0 data. The only difference from Fukagawa et al. 
(2013) was the visibility weighting applied in the final step of the imaging; we adopted in this study the Briggs weighting with a 
robust parameter of 0.5 to best recover the weak and extended components of the emission. The size in FWHM and the position angle 
for the major axis of each synthesized beam were 0('47 x 0!'40 = 60 AU x 56AU {PA — 59.9°), 0!'50 x 0"42 = 70AU x 59AU 
(PA = 57.4°) and 0"50 x 0"42 = 70AU x 59AU (PA — 60.6°), for the continuum at 336 GHz, ^®CO and C^®0, respectively. 
Further details on bandpass and gain calibrations are described in Fukagawa et al. (2013). The rms noise is 0.13 mjy beam“^ for 
the continuum whereas it is 6.1 and 8.3 mJy beam“^ in the 0.12 km s“^ wide channels for the line emission of ^^CO and C^®0, 
respectively. 


2.2 Results of Dust Continuum Emission 


Figure 1 shows the continuum map, which looks quite similar to that presented by Fukagawa et al. (2013). The outer asymmetric 
ring as well as an inner unresolved component are detected, and these are separated by a radial gap. The position of the unresolved 
component coincides with the velocity centroid of ^^CO J = 3 — 2 (section 2.3), and it is regarded as the stellar position in the 
following. The radial profiles of surface brightness is well described by a Gaussian function, and the brightest and faintest of their 
peaks are located at PA « 23° and PA « 223°, respectively (see also figure 4 of Fukagawa et al. (2013)). The averaged surface 
brightness profiles I{r) in PA = 11° — 31° and 211° — 231° are fitted by 


I{r) = Ip exp 


iX t‘0,obs) 


( 1 ) 


where r is the angular distance from the star, ro.obs is the peak position and Wohs is the width of the Gaussian. The best-fit profiles 
are shown in figure 2 and the best-fit parameters are summarized in table 1. As shown in table 1, the contrast in Ip is 24 between 
these two position angles. In the following, we refer to the averaged profiles in fhe section of 11° < PA < 31° as “north profile” 
and that in 211° < PA < 231° as “south profile” (see the right panel of figure 1). 


2.3 Results of and C^®0 J = 3 - 2 Emission 

Figures 3 and 4 show the moment maps of the J = 3 — 2 and C^®0 J — 3 — 2 emission, respectively. Moments 0, 1, and 
2 correspond to the integrated intensity, intensity-weighted mean velocity, and velocity dispersion, respectively. As shown in the 
appendix, the emission above 5 ct level is detected in ulsr = (0.64 — 7.00) km s“^ in and wlsr = (1.24 — 6.16) km s“^ in 
C^®0. The azimuthal asymmetry is weak in the moment 0 map although the northern part tends to be slightly weaker, possibly due 
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Fig. 1. (Left) Map of the continuum emission at 336 GHz with the Briggs weighting with the robust parameter of 0.5. The synthesized beam, Of'47 x 0('40 
with the major axis PA = 59.9°, is indicated by the white eliipse in the bottom left corner. The contours correspond to 1, 10 and 100 mjy beam“^. The Icr 
levei is 0.13 mJy beam^^ (Right) The position angie of the major axis is indicated by dashed iine. The regions where the azimuthal average is taken to obtain 
the north profiles {PA = 11° - 31°) and the south profiles {PA = 211° - 231°) are indicated by red and blue hatches, respectively. Contours are the same 
with the left panel. 


r[AU] 

0 100 200 300 400 



r[AU] 

Fig. 2. Radial profiles of surface brightness of the continuum emission at 336 GHz in PA = 11 - 31° (top panel, squares) and PA = 211 - 231° (bottom 
panel, squares). Solid lines indicate the best-fit Gaussian function (see Table 1 for parameters). Error bars indicate the standard deviation after the averaging 
over 20 degrees in PA. 
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Fig. 3. Moment maps of the J = 3 - 2 line, (left) Moment 0 map, Integrated over the velocity range of ulsr = (0.40 - 7.00) km s ^. The white contour 
shows the 5cr ievel (48 mjy beam'^ km s~^). (middle) Moment 1 map, created by the emission above the So- levei in channel maps of 0.12 km s~^ resolution 
(figure 32-34 in appendix). The contours along PA « 71° are those of the systemic velocity (dlsr = 3.7 km s^^^), and the contour spacing is 0.5 km s“^. 
(right) Moment 2 map, created by the emission above the So- level in channel maps of 0.12 km s'^ resolution. Contour starts at 0.2 km s“^, and its spacing 
is 0.2 km s~^. The synthesized beam, o."50 x o."42 with the major axis PA = 57.4°, is indicated by the ellipse in the bottom left corner of each panel. 


[Jy beam"' km s"'] [km s'^] dV [km s"'] 

-0.05 0 0.05 0.1 0.15 0.2 0.25 2 3 4 5 0 0.1 0.2 0.3 0.4 0.5 



dR.A. [arcsec] 


Fig. 4. Moment maps of the J = 3 - 2 line, (left) Moment 0 map, integrated over the velocity range of dlsr = (l.O - 6.40) km s The white contour 
shows the 5<t level (55 mJy beam“^ km s“^). (middle) Moment 1 map, created by the emission above the 5tT level in channel maps of 0.12 km s“^ resolution 
(figure 35-36 in appendix). The contours along PA ra 71° are those of the systemic velocity (dlsr = 3.7 km s”'^), and the contour spacing is 0.5 km s“^. 
(right) Moment 2 map, created by the emission above the 5<t level in channel maps of 0.12 km s~^ resolution. Contours at 0.2 and 0.4 km s~^ are shown. 
The synthesized beam, o."50 x o."42 with the major axis PA = 60.6°, is indicated by the ellipse in the bottom left corner of each panel. 


to higtier continuum levels. Furthermore, the line emission (^^CO, in particular) is clearly detected in the inner regions down to 
r « 20 AU (0!'15) as well as the outer regions out to r « 400 AU (2"8). The velocity distribution revealed in moments 1 and 2 
is consistent with a disk in Keplerian rotation. A constant (systemic) velocity of 3.7 km s“^ is found along PA = 71° ± 2° in the 
moment 1 maps, and this is regarded as the direction of the minor axis of the system. The position-velocity (P-V) diagram along 
the major axis {PA = —19°) is explained well by Keplerian rotation with stellar mass of 2.2Mq and the inclination angle of 27°, 
as described in detail in the appendix (see also Fukagawa et al. 2013; Perez et al. 2015). We adopt these values for PA of the major 
axis and inclination of the system throughout this paper. 

Figure 5 shows the north and south profiles of the moments 0 and 1, which will be the main focus of the modeling described in 
later sections. It is clear that the moment 0 profiles of both and €*^*0 are very different from those of dust continuum emission; 
these are more extended than the Gaussian-like dust emission profiles. The moment 1 profiles of these two lines agree with each 
other in lOOAU <r< 280AU, indicating that both these lines successfully reveal the Keplerian rotation in the regions where the 
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Fig. 5. The north (top row) and the south (bottom row) radial profiles of the observed moment 0 (left panels) and moment 1 (right panels) of (red circles) 
and (blue squares). The systemic velocity of 3.7 km/s is subtracted in calculating the moment 1 profiies. Error bars indicate the standard deviation after 
the averaging over 20 degrees in PA. The radial profiles of the continuum emission is overplotted in the moment 0 radial profiles for comparison. 


emission is detected with a high signal-to-noise ratio (S/N). The moment 1 profile of further reveals the gas motion down to 
r « 20 AU and out to r « 400 AU. As discussed in the subsequent sections, this indicates the existence of tenuous gas components in 
these inner and outer regions. The line has larger moments 0 and 2 in almost all the positions than the C^®0. This is because 
has a larger optical depth at every velocity channel and hence has a broader line profile than C^®0 (see also section 6.2.2) 


3 Method of Modeling 

Our goal is to find densify and temperafure distribution of the disk around HD 142527 based on the continuum and CO line emission 
profiles described in the previous section. We construct axisymmetric models that can reproduce similar radial brightness profiles in 
the direction where dust emission is the brightest (north profile) and the faintest (south profile) since the model that fully accounts 
for the azimuthal asymmetry can be very complex. This approach is similar to that taken by Bruderer et al. (2013) in making the 
model for the disk around Oph IRS 48, which also exhibits strong azimuthal asymmetry. 

3.1 Dust Distribution Modeis 

We first derive dust density and temperature distributions from the continuum emission under the assumption that the disk is in 
thermal and hydrostatic equilibrium. For simplicity, the gas-to-dust ratio is assumed to be uniform in the vertical direction and the 
gas temperature is assumed to be the same as dust. In other words, we have ignored dust sedimentation (e.g., Dubrulle et al. 1995) 
and temperature difference between gas and dust in the low density upper layers (e.g., Kamp and Dullemond 2004; Nomura and 
Miller 2005). We note that the gas-to-dust ratio can vary in the radial direction, which is the main focus of the modeling of gas 
observations (see sections 3.2 and 5). The disk is assumed to be heated only by stellar irradiation, since viscous heating is less 
important in the region considered in this paper. The star is assumed to have the effective temperature of Ten = 6250 K and the 
radius of i? = 3.8 Rq (Verhoeff et al. 2011).* We solved the zeroth and first order moment equations of the radiative transfer (Ml 
method, see Kanno et al. 2013). We use 226 colors in the wavelength range of 0.1 pm < A < 3.16 mm, resulting in the spectral 
resolution of AlogA = 0.02 (i.e., A/AA = 21.7). The computational box covers 30 AU <r< 410 AU and \z\ < 120 AU with the 
spatial resolution of 2 AU in the cylindrical coordinates. The disk is assumed to be symmetric with respect to the midplane. 

The dust is assumed to consist of silicate, carbonaceous grains, and water ice having the mass fractional abundance of (^su = 
0.0043, Ccarbon = 0.0030, and (^ice = 0.0094, respectively, which are consistent with solar elemental abundance (Anders and Grevesse 

' The stellar parameters are updated in Mendigutia et at (2014), but the stellar luminosity is within the error for the parameters described here. 
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Fig. 6. The absorption (red) and the effective scattering (blue; see text for definifion) coefficienfs for the dust model used in this paper. 


1989). The dust particles are assumed to have the power-law size distribution of oc ® with the maximum size of Umax = 1 mm 
(Nomura and Miller 2005). Figure 6 shows the absorption (Ki.,a) and the effective scattering coefficients in units of cm^ per unit 
gram of dust. Here, we define the effective scattering coefficient as (1 — g), where is the scattering coefficient and g = (cos9) 
is the scattering asymmetry factor.^ The model with the maximum dust size of 1 mm gives the large value of opacity at sub-mm 
range (Aikawa and Nomura 2006). Consequently, the dust mass evaluated in this paper should be close to the lowest. Our opacity 
at 870 nm has values of Ka = 2.9 cm^/g with /3 ~ 1.2 at 0.3 mm < A < 1 cm and «:s(l — g) = 26 cm^/g per unit dust mass. The 
absorption coefficient at 870 fim is 20% smaller than the value adopted by Beckwith et al. (1990), which is ~ 3.5 cm^/g per unit 
dust mass. We also note that the effective scattering coefficient is large compared to the absorption one in the sub-mm range since 
the maximum dust size is comparable with the wavelength. 

The dust surface density is assumed to have the form 


Ed(r) = Ed,0 exp 



( 2 ) 


since the radial profile of the continuum emission is well approximated by a Gaussian function. For an assumed set of parameters 
(Sd,o, Td, uid), we obtain the spatial distribution of dust density pd{r,z), temperature T{r,z), and the radiation energy density 
Jv{r, z) for the 226 colors at each grid cell. 

The expected observed surface brightness profiles of dust continuum emission is obtained by ray-tracing. The surface brightness 
7^ at frequency iz is calculated by solving 

^ = -PdX. [/. - S.], (3) 

as 

where s is the coordinate along the line of sight, Xv is the extinction by absorption and (effective) scattering, Xi^ = + (1 — g)>^v,s. 

The source function, S^, is given by 


Sv ~ iX — uJiiJv, (4) 

where = (1 — g)K^,slXv is the (effective) albedo^ and Bv{T) is the Planck function. The second term in equation (4) represents 
the scattered light, and the scattering is assumed to be isotropic for simplicity. The optical depth Ti, along the line of sight is given 
by 

^ = -X.P.- (5) 

We compute the model image by using {Nii,Ng,) = (128,128) rays covering 35 AU < R < 400 AU and 0 < 'F < 27r region, where 
(7?,T>) are the polar coordinates on the sky-plane with the central star at the origin. We use the inclination angle of 27°, as described 
in section 2 


- We have left the term (1 - g) for consistency with the Ml method. See Mihalas and Mihalas (1984); Gonzalez, Audit and Huynh (2007) for the appearance 
of (1 — 9 ) factor in the Ml method. 

^ With the dust model of o„ax = 1 mm, the value of Qu = K.s/{na. + k.^) and given in the main text differs only by ~ 5% at the wavelengths of interest. 
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For comparison with observations, we convolve the model images with the Gaussian function of the same beam size and ori¬ 
entation with the observations. We extract the radial surface brightness profiles from the convolved image and compare them with 
observations. We iterate this procedure until the given set of parameters {T,d,o,rd,Wd) reproduces the observed profiles shown in 
figure 2 reasonably well. 


3.2 Gas Distribution Modeis 


We then derive the gas distribution that accounts for both moments 0 (integrated intensity) and 1 (intensity-weighted mean velocity) 
radial profiles in the north and south directions. We have chosen these moments of the line emission because they are least affected 
by beam dilution. It should also be noted that the observed moment 2 profiles contains the uncertainty coming from the choice of 
cutoff levels when producing the moment map from the data. 

In later sections, we show that uniform gas-to-dust ratio models do not reproduce the observations. Therefore, we assume that 
the gas density pg (r, z) is given by 

Pg{r,z)=£,{r)pd{r,z), ( 6 ) 

where ^(r) represents the gas-to-dust ratio at each radius. The gas surface density 'Sg{r) is then given by 

E3(r)=C(r)E,(r) (7) 


In later sections, we look for the forms of ^(r), or, equivalently, the form of Eg(r), that best reproduces the observed radial profiles 
of momenls 0 and 1. 

We assume that the gas rotation is Keplerian at the disk midplane, 

«.ot(r) = 3.13(^^) ' [kms-i], (8) 

where 3.13 km s“^ Is the Kepler velocity at 200 AU around a 2.2 Mq star. Rotation velocity can be slightly different from Keplerian 
due to, for example, radial pressure gradient force, but the difference is at most of the order of thermal velocity (~ 10% of the Kepler 
rotation velocity), which is hardly observed with current velocity resolution. Gas temperature T{r,z) is assumed to be the same as 
the dust temperature as mentioned in section 3.1. 

The expected brightness of line emission is calculated with ray-tracing methods by solving 

^ = -XuiiIui-S^i), (9) 

as 

where lui is the intensity of the line emission from the upper state u to the lower state L The total extinction Xui comes from both 
dust and gas, 


\ul PdXiy T ^ul . t 

47r 

and the source function Sui is given by spontaneous emission and dust emission Si^ given in equation (4), 

£? _ 1 _ ^ ^ I pdXi' c- 

^ul - ^u-^ul^ul . \ Ju' 

Xui 47r Xui 


( 10 ) 

( 11 ) 


Here, n„ and ni are level populations for the upper and lower state, respectively, Bui, and Biu are Einstein coefficients, and 
4>ui is the line profile function. We note that the scattering of line emission by dust particles is not included in this work. Local 
thermal equilibrium (LTE) is assumed to calculate the level population, which is a valid assumption for lower transition lines of CO 
in a protoplanetary disk where typical density is much higher than the critical densities for these lines (Pavlyuchenkov et al. 2007). 
The fractional abundance of and C^®0 is assumed to be 9 x 10“^^ and 1.35 x 10“^ (Qi et al. 2011). We have assumed that the 
line width is determined by thermal broadening. 

We construct the model channel maps from —4.5 km s“^ to 4.5 km s“^ with respect to the systemic velocity with 0.06 km s“^ 
step and each model channel map is convolved with the Gaussian beam. The model radial profiles of the moments 0 and 1 are then 
extracted in the north and south directions for comparison with the observations. The continuum emission is subtracted on the image 
by using line-free channels and the velocity channels showing emission only below the detection limit are excluded when calculating 
the moment maps from the model. To confirm the validity of continuum subtraction in the image plane, we have also made imaging 
simulations for several cases in a more rigorous way in which the visibilities for the baselines sampled in our observations are first 
generated from the disk model with the CASA simulator and then continuum subtraction is made in the uu-plane. After imaging 
with the same procedure as that described in section 2.1, we have confirmed that the resultant line profiles agree wilh those obtained 
with the image-based continuum subtraction within a few percent. 
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Fig. 7. The comparison of continuum emission and the modei surface brightness for the north (top) and south (bottom) profiles. 


4 Results of Dust Distribution 

We have searched for the parameters Ed,o> Td, and Wd in equation (2) that can reproduce the observed north and south profiles. 
In the south profiles, the parameter search is straightforward since the dust emission is optically thin. In the north, on the other 
hand, it is necessary to carefully look at the dependence of surface brightness profiles on dust distribution parameters since the dust 
emission is optically thick (see section 6.2.1). The details of the parameter search for the north profiles are summarized in appendix 
2. The best-fit results are summarized in table 2 and figure 7 compares the radial profiles of the dust continuum emissions derived 
from the model and observations. The best-fit models have the dust density and temperature distributions shown in figure 8. Our 
derived parameters for the dust disk (or ring) radius and width is consistent with the results of Verhoeff et al. (2011), who show that 
the massive outer disk extends from 130 AU to 200 AU from the central star based on their modeling of the SED and mid-infrared 
images (and therefore dust emission). 

The dust ring emission is only marginally resolved in ALMA Band 7 since the full width at half maximum (FWHM) of the 
radial Gaussian function of the model surface density is ~ 2y/ln2wd ^ 50 AU (see table 2), which is slightly smaller than the beam 
size (~ 0"45 ~ 60 AU at 140 pc). As a result, the radial width of the observed surface brightness profile (~ 2'\/ln2M;obs ~ 85 AU 
in FWHM; see table 1) is larger than that of the model surface density. Perez et al. (2015) measured the radial width of the dust 
continuum emission at 230 GHz to be 0"9, while at 345 GHz to be ~ 0!'55. The difference of the width at different band may be 
explained by the effect of the convolution by the beam. The radial width of the continuum emission is not well resolved. The ratio 
of the radial width of the continuum emission between 230 GHz and 345 GHz is ~ 1.6, which is close to the ratio of the frequency 
(and thus the ratio of the beam size) between the two bands. 

It is indicated that the amount of dust particles is ~ 70 times more at the peak of the north profiles than at the peak of the south 
profile, although the surface brightness contrast of the dust emission between the north and the south peak locations is ~ 24. The 
difference between the surface density contrast and the surface brightness contrast is largely due to the fact that the northern part of 
the disk is optically thick to dust emission. It should also be noted that the scattered light component in the dust continuum emission 
is not negligible in the north profile (see section 6.2 for discussion). 
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Fig. 8. Dust density and temperature distribution for the best-fit models of north (top) and south (bottom) profiles. Gas density is indicated by color in unit of 
g/cm^ and the temperature is shown by contours 



North 

South 

Ed,o [g cm 

0.6 

8.45 X 10"® 

rd [AU] 

173 

196 

Wd [AU] 

27 

34 


Table 2. Best-fit values for dust distribution. 

5 Results of Gas Distribution 

We now turn our attention to gas distribution. We first show that the models with uniform gas-to-dust ratio 100 do not reproduce 
the observed profiles in section 5.1. We then describe in detail how we construct the models for gas distribution step by step in 
subsequent subsections. Our final results are summarized in section 5.6. 

5.1 Failure of Uniform Gas-to-Dust Ratio Models 

We first show that the models with uniform gas-to-dust ratio 100, i.e., ^(r) = 100 in equation (6), fail to reproduce the observed 
radial profiles of momenfs 0 and 1. In this case, the radial gas surface density profile Eg(r) = ^(r)Ed(r) is given by a Gaussian 
function as the dust distribution is. Figure 9 shows the radial profiles of moments 0 and 1 for the north direction and figure 10 shows 
the same but for the south direction. 

At inner radii (r < 100 AU), the model profiles of moment 0 show too weak emission compared to observations and that of 
moment 1 tends to be too slow compared to the observations. These indicate that there has to be more gas within 100 AU. On 
the other hand, around the peak of the dust distribution (r ~ 170 AU for north and r ~ 200 AU for south), the model moment 0 
is too bright compared to the observations, especially in the case of €*^*0 profiles. This suggests that gas-to-dust ratio is smaller 
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Fig. 9. The radial profiles of moment 0 (left panels) and moment 1 (right panels) for (top row) and (bottom row) for the north profiles in the case 
where gas-to-dust ratio is fixed to 100 everywhere in the disk. 
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Fig. 10. The same as figure 9 but for the south profiles. 
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Eo [g/cm^] 0.0845, 0.2325, 0.845, 2.325 

To [AU] 90,100,110,120 

rout[AU] 200,250,300,350 

Table 3. Parameters explored for power-law gas distribution given by equation (12). 

than 100 at places where dust particles are concentrated. It is impossible, however, to decrease the amount of gas at ~ 200 AU 
and to increase at inner radii simultaneously under the assumption of uniform gas-to-dust ratio. The dust distribution is already 
determined in section 4. Therefore, we need to consider completely different radial distribution of the gas from the dust by varying 
the gas-to-dust ratio ^(r) within the disk. 

5.2 Power-Law Gas Profiles 

Having found that the uniform gas-to-dust ratio models do not reproduce the observed profiles of gas emission, we now explore 
models of gas distribution, which is not necessarily restricted to constant gas-to-dust ratio. To acquire consistent results with the 
dust distribution calculations described in section 3.1, we still assume that the gas-to-dust ratio is constant in the vertical direction at 
each radius, but it is not constant in the radial direction. With a trial function of gas surface density ^g{r), it is possible to calculate 
the gas-to-dust ratio at each radius by calculating ^(r) = Eg(r)/Ed(r), which is then used to obtain the gas density at each grid cell 
from equation (6). 

From the results presented in section 5.1, it is indicated that the gas distribution is more extended in the radial direction than dust. 
Therefore, we first try the model in which the gas surface density profile Eg(r) is given by the r~^ power-law profile at rc < r < rout, 

(^c<r<rout), (12) 

and zero otherwise. Here, can be interpreted as the cavity radius. Since constant gas-to-dust ratio models fail to explain the 
observed profiles, it is necessary to try models with different gas distribution. The gas-to-dust ratio is not constant in such cases. 
We have first tried several models with Gaussian surface density profiles with different parameters, but we have found that these 
models fail to reproduce the observed profiles. This is primarily due to the fact that the gas emission extends to large radii, but 
Gaussian function falls off too rapidly. Therefore, we consider power-law distribution of the gas surface density, which is more 
radially extended than Gaussian profiles. We choose power-law distribution as a representative model of such radially extended 
gas profiles and fix the power-law exponent to —1 in order to reduce the number of free parameters. The choice of r~^ power is 
motivated by the models of steady state accretion disk, but we do not claim that this is the only possible solution. We aim to find 
one possible disk model that reasonably reproduces observations and discuss its implications for the overall picture of how the dust 
and gas are distributed. 

We explore the parameter space of (Eo,rc,rout) and look for values that match the observed moments 0 and 1 profiles for the 
north and south profiles. Our strategy is to find first the values of Eo that reasonably matches with the observations at 150 AU < 
r < 200 AU, where dust emission is bright. Then we look for the values of by investigating the inner radii. Finally, we search for 
the values of Tout by investigating the outer radii. We show the results one by one in subsequent sections. The values of these three 
parameters explored in our set of calculations are summarized in table 3. 

5.3 The Overall Gas Distribution 

We first find the values of Eo that reasonably match the observed profiles at 150 AU ^ r ^ 250AU. For this purpose, we fix 
(fc, font) = (100 AU, 280 AU) for both north and south profiles and we vary Eo from 0.0845 g cm“^ to 2.325 g cm“^ with a step 
of a factor of ~ 3 (see table 3 for specific parameters). The lowest value of Eo is chosen in such a way that the gas-to-dust ratio 
at the peak of the dust distribution in the south direction is 10. Figures 11 and 12 show the results for the north and south profiles, 
respectively. 

For the north profiles of C^®0 moment 0, the observations fall between the models with Eo = 0.845 g cm“^ and Eq = 
2.325 g cm“^. For moment 0, the model with Eq = 2.325 g cm“^ seems to show better fit, but it is actually difficult to 
judge which is plausible when the existence of remaining gas inside the cavity is considered (see section 5.4). Note that all the 
models fail to explain the moment 0 profile within the cavity (r < rc = 100 AU) as the model profiles fall off too rapidly towards 
inner radii. This is because there is no gas emission inside 100 AU in the model. We also point out that the model moment 1 profiles 
show too slow velocity at inner radii. This is because moment 1 is calculated to be zero in such gas-free regions before we convolve 
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Fig. 11. The radial profiles of moment 0 (left panels) and moment 1 (right panels) for (top row) and (bottom row) for the north direction when the 
parameter So is varied. 
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Fig. 12. The same as figure 11 but for the south profiles. 
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the model with the beam. The apparent emission inside the cavity in the model is the result of beam dilution. Similar discrepancy 
between the model and observations can be found at the outer radii of r ^ 300 AU, while in this case the beam dilution is not as 
significant as inner radii. 

One prominent feature of the model north profile is the two peaks of the surface brightness profiles at r ~ 100 AU and r ~ 250 AU, 
while the gas surface density profile is smooth. These peaks are primarily due to the existence of large amount of dust in relatively 
narrow radial range. The dust emission is optically thick and is very bright. The line emission is partially hidden by the optically 
thick dust. Therefore, the difference in brightness between dust continuum and line emissions is not very significant. The line 
emission at the peak of the dust emission is largely affected when the continuum emission is subtracted to calculate moment maps. 
The two “peaks” of gas emission should be considered as “one trough” caused by the subtraction of bright continuum emission. 
Although not as prominent as the model profiles, it is possible to see this effect in the observed profiles, which show a slight dip of 
the moment 0 at ~ 170 AU (see figure 5 and 11). We further discuss about these apparent bump and/or trough structures in section 
6.3. 

For the south profiles, Eq = 0.2325 g cm“^ models show reasonable fit to both and C^®0. However, the discrepancy 

between the model moment 1 profiles and the observations at inner and outer radii is present, as in the north profile case. This 
discrepancy is further studied in sections 5.4 and 5.5. We note that the “trough” structure is not very significant in the south profile 
since the continuum emission is weak and the continuum level is sufficiently low. Comparing the values of Eq derived for the south 
profile with that for the north, it is indicated that the azimuthal asymmetry of gas distribution is not as significant as dust distribution. 
The azimuthal contrast of gas surface density is only by a factor of ~ 3 — 10 between the north and south profiles, while the peak 
dust surface density differs by a factor of ~ 70. 

5.4 Inner Radius of the Gas and Remaining Gas inside Cavity 

We now explore the parameter space forrc, the cavity radius of the gas. We fix (Eo,rout) = (0.845 g cm“^,280 AU) and (Eo,rout) = 
(0.2325 g cm“^, 280 AU) for the north and south profiles, respectively. We vary rc from 90 AU to 120 AU, and see whether we 
can match the observed profiles within r ^ 150 AU. Especially, we look for values of which matches the peak locations of the 
moment 0 profiles residing at r ~ 110 AU in the north profile and at r ~ 140 AU in the south profile. 

Figures 13 and 14 show the results for the north and south profiles, respectively. It is shown that Vc ~ 100 — 110 AU can 
explain the location of the peak in the profiles of moment 0. However, all the models give slower velocity in moment 1 profile than 
observations at r < 100 AU. It should also be noted that the model profiles of moment 0 fall to zero quickly towards the inner radii 
while the observed profiles, especially south profiles, show more gradual decrease. 

The discrepancy between the model and observed profiles at inner radii can be explained if we assume that there is some re¬ 
maining gas at r < rc. The existence of gas inside the cavity is also indicated by the existence of emission all the way to 
the central star (Casassus et al. 2013; Rosenfeld et al. 2014; Perez et al. 2015). To explore the parameters, we assume that the gas 
surface density at r < rc is given by 

E,(r) = /i„E,(rc)(^) (r < rc) (13) 

while the gas at r > rc is kept the same as equation (12). The parameter ./■in is varied from 0 (complete cavity) to 3/16 for the north 
models and from 0 to 3/8 for the south models (see table 4). Figures 15 and 16 show the results for the north and south profiles, 
respectively. It is clear that the match between the model and observations is better for the models with remaining inner gas. The 
values of /in that give the best match between the model and observations seem to be ^ 1 /8 for the south profile. For the north 
model, the best-fit lies somewhere between 1/16 and 1/8. To keep the parameter search simple, we allow a factor of ~ 2 error here, 
and use 1/8 as a representative value. An important indication from the modeling approach is that similar gas distribution models 
(within a factor of 2 — 3) can account for the observations both in the north and south regions of the disk, which is very different 
from the case of dust distribution. 

We have assumed that the gas surface density increases linearly with radius within r < rc. However, the functional form of the 
gas distribution does not affect the results. We have checked this by using a series of models with constant gas surface density at 
r < Tc. This is partly because the beam size is relatively large and the beam dilution effect is significant especially when we discuss 
the gas distribution at inner radii. We therefore conclude that it is a robust conclusion that there should be some remaining gas 
within the cavity but there is an uncertainty in the details of how the gas is distributed. Indeed, Perez et al. (2015) have used different 
functional form for the models of gas inside the cavity and obtained results consistent with observations. The amount of gas within 
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Fig. 13. The radial profiles of moment 0 (left panels) and moment 1 (right panels) for (top row) and (bottom row) for the north direction when the 
parameter is varied. 
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Fig. 14. The same as figure 13 but for the south profiles. 
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Eo [g cm^] (north) 0.845 

rc [AU] (north) 110 

rout[AU] 280 

/in (north) 0, 1/16, 1/8, 3/16 


Eo [g/cm^] (south) 0.2325 

rc [AU] (south) 100 

rout[AU] 280 

/in (south) 0, 1/8, 1/4, 3/8 


Table 4. Parameters explored for remaining gas inside cavity 
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Fig. 15. The radial profiles of moment 0 (left panels) and moment 1 (right panels) for (top row) and (bottom row) for the north direction when the 
parameter /i„ is varied. 


the cavity in Perez et al. (2015)^ is ~ 1 — 2 x 10“® Mq while we have the gas mass inside the cavity being ~ 5 x 10“^ Mq if we 
integrate the model with (Eo,rc,/in) = (0.845g/cm^, 110 AU, 1/8), which is preferred for the north profile, from r = 0 to r = rc. 
These results are consistent within an order of magnitude. We consider that the difference comes from the fact that they use lower 
temperature and smaller line width than our model in estimating the gas mass inside the cavity (see also section 6.2.3). 


5.5 Outer Radius of Dense Gas and Remaining Tenuous Gas at Large Radii 

We now turn our attention to the gas distribution at outer radii. In the series of models presented in sections 5.3 and 5.4, the outer 
radius of the gas distribution is fixed to rout = 280 AU. As indicated in figures 11-14, the moment 1 radial profiles of cuts off 
at ~ 350 AU while the observations indicate that there must be some gas at outer radii. 

To explore the parameter space, we first vary rout while keeping (Eo,rc,/in) = (0.845 g/cm^, 110 AU,0) and (Eo,rc,/in) = 
(0.2325 g/cm^, 100 AU, 0) for the north and south models, respectively. We have checked that the details of the choice of rc and /in 
does not affect the brightness profiles at outer radii. Figures 17 and 18 show the results for the north and south profiles, respectively. 
It is shown that the moment I profiles indeed extend to larger radii as we increase rout. However, if rout exceeds 250 AU, the models 
show much larger values of moment 0 than observed at outer radii. 

The discrepancy of the moment 0 profiles between the observations and models indicates that we need to decrease the amount of 
gas at outer radii than currently assumed power-law profiles, especially at r > 250 AU. At the same time, there should still be some 
gas remaining at outer radii so that moment 1 profiles do not cut off. As an alternative series of models, we now try commonly used 
exponential taper model given by 

^•M = ^«(20i(Au)"“-(-/T) 


The cavity radius is 90 AU in their modei 
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Fig. 16. The same as ligure 15 but tor the south protiles. 
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Fig. 17. The radial protiles ot moment 0 (left panels) and moment 1 (right panels) lor (top row) and (bottom row) tor the north direction when the 
parameter rout is varied. For this parameter search, we assume that there is no gas outside rout. 
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Fig. 18. The same as figure 17 but for the south profiles. 
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Fig. 19. The radial profiles of moment 0 (left panels) and moment 1 (right panels) for (top row) and (bottom row) for the north direction when the 
parameter rout is varied. For this parameter search, we assume that the gas at outer radii exhibits exponential cutoff as in equation (14). 


for r > Tc. The results are shown in figures 19 and 20 for the north and south profiles, respectively. Although the model moment 
1 profiles show good match with observations for both the north and south profiles, the models still exhibit too large values for the 
moment 0 profiles. This indicates that the gas surface density should decrease more rapidly than exponential profiles at outer radii. 

To incorporate the rapid decrease of gas surface density at outer radii, we finally try a series of model gas surface density profiles 
fhat cuf off at rout = 250 AU but there is remaining, tenuous gas with constant surface density at r > rout all the way to the outer 
edge of the computational domain. In this framework, the gas surface density at r > rout is given by 


^g,out (r) - /out Sg('^out) (j" ^ ^^out) 


(15) 
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Fig. 20. The same as figure 19 but for the south profiles. 


Eo [g/cm^] (north) 

0.845 

Eo [g/cm^] (south) 

0.2325 

rc [AU] (north) 

no 

rc [AU] (south) 

100 

rout[AU] 

250 

'Tout [AU] 
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/out 

1 , 10-\ 10■^ 10"® 

/out 

1 , 10-\ 10■^ 10”^ 


Table 5. Parameters explored for tenuous gas outside 


where Eg (rout) is the surface density at r = rout given by the power-law distribution in equation (12) and /out controls the amount 
of remaining gas at outer radii. The parameters explored are given in table 5. 

Figures 21 and 22 show the results for the north and south profiles, respectively. For the north profiles, the model profiles 
at r > 250 AU show reasonable match with observations when /out ~ 10“^. For the south profiles, /out = 10“^ models show 
agreement at r < 300 AU while they show brighter emission than observations at the outermost radii. The moment 0 profiles of the 
model with /out = 10~^ is fainter than observations all the way in the outer radii. Therefore, the actual values of /out may vary as a 
function of radius for the south profiles, but the overall value of /out may be around ~ 10“^ to ~ 10“^. Finally, we briefly note that 
the appearance of the “trough” at r ~ 170 AU in the north profile (see discussions in section 5.3) is different when we use different 
profiles at r > Tout. This is because the change in gas surface density at outer radii (within the parameters explored in this section) 
mainly affects the moment 0 profiles at r > 200 AU. The strength of the second (apparent) bump at r ~ 250 AU is therefore affected. 


5.6 Summary of Gas Distribution 

We have looked for the gas distribution models that reproduce the observed moments 0 and 1 profiles of and C^®0. We 
have shown that the gas-to-dust ratio is not constant in the disk and that the gas distribution should be more extended in the radial 
directions than the dust distribution. We have seen that there should be some remaining gas inside ~ 100 AU and some tenuous gas 
at r ^ 250 AU. 

To be more specific, we have assumed that the model surface density of gas is given by 

/incr(rc) (r < rc) 

o-(r) (rc < r < rout) ’ 

/out '^('^out) ^ ^out) 

where a(r) is given by the power-law distribution, 
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Fig. 21. The radial profiles of moment 0 (left panels) and moment 1 (right panels) for (fop row) and (bottom row) for fhe north direction when the 
parameter /out Is varied. For this parameter search, we assume that the gas at outer radii exhibits a cutoff at 250 AU, but there is still some remaining gas 
outside the cutoff radius as in equation (15). 
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Fig. 22. The same as figure 21 but for the south profiles. 
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Table 6. Parameters of the reference models for gas distribution. 
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Fig. 23. North radial profiles of moment 0 (left panels), moment 1 (middle panels), and moment 2 (right panels) of (top row) and (bottom row) for 
the reference model given by equation (16). The values of the model parameters are given in table 6. 


There are five control parameters: Eo that determines the overall amount gas surface density, Tc that determines the inner cutoff 
(cavity) radius, /in that determines the amount of gas within the cavity. Tout that determines the outer radius of the disk, and /out 
that determines the amount of tenuous gas at large radii. The values of these parameters that give reasonable fit to observations are 
given in table 6. Hereafter, we call these best-fit models “reference models”. 

The comparison of the radial profiles between the models and observations is given in figures 23 and 24 for the north and south 
directions, respectively. For these models, we show the comparison between the models and observations of the moment 2 radial 
profiles averaged over 11° < PA < 33° (north) and 211° < PA < 233° (south) as well as the moments 0 and 1, which have been 
the main focus of the modeling. The radial profiles of all the moments show reasonable match between observations and models. 
Finally, we show the comparison of P-V diagrams in the north and south directions for the reference models in figure 25. The models 
and observations show reasonable match not only in the radial profiles of the moment maps but also in the P-V diagram. 

The radial profiles of the model moment 0 profiles exhibit some bumps and troughs while observations show more smooth 
profiles. This is most prominent at around ~ 170 — 200 AU region of the north profiles where the dust emission is very bright. These 
bumps are apparent structures due to the subtraction of bright continuum emission, and we discuss about this further in section 6.3. 


6 Discussions on the Reference Model 

We have constructed the models of dust and gas distribution in the protoplanetary disk around HD 142527. In this section, we look 
at our reference models of dust and gas distribution and discuss indications and caveats of our model. 
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Fig. 24. The same as figure 23 but for the south profiles. 
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Fig. 25. Comparison of P-V diagram between observations (black contours) and the model (red contours) for the north and south directions. The left panels 
show the north direction and the right panels show south direction, while the top row shows the results for and the bottom row shows Contours 
start from 5-<t level (0.135 Jy/asec^ for '^^CO and 0.175 Jy/asec^ for C^®0, respectively) with lO-o- step. 
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Fig. 26. Surface density distribution of gas (solid iine) and dust (dashed iine) in northern (red) and southern (blue) part of the disk. The left panels show 
surface density in the linear scale (top panel for the north and the bottom panel for the south model), while the right panel shows the surface density in the log 
scale. 

6.1 Summary of Gas and Dust Distribution 

We have shown that the gas and dust are distributed very differently in the disk. Figure 26 shows the surface density profiles of 
gas and dust in the north and the south regions for the reference model. Dust distribution can be explained by the radial Gaussian 
ring-like profile with the width of Wd ~ 30 AU (50 AU in FWHM) for both the northern and southern regions. Flowever, there is 
a factor of ~ 70 difference in surface density between the two regions (see equation (2) and parameters listed in table 2 for dust 
distribution). Gas distribution is more or less axisymmetric and radially extended. The variation of surface density between the 
northern and southern part is a factor of ~ 3 and the gas density cuts off at ~ 100 AU, and Tout ~ 250 AU, but both inner and 
outer regions are not completely devoid of gas (see equation (16) and parameters listed in table 6 for gas distribution). 

The gas-to-dust ratio is ~ 30 at the peak of the dust emission in the southern part of the disk, while it reaches ~ 3 in the northern 
peak. Since dust emission is optically thick in the northern part of the disk, the dust emission is not very sensitive to the amount 
of dust (see appendix 2) so the amount of dust in the northern part may be considered as a lower limit. Therefore, the gas-to-dust 
ratio may be even lower in the northern part of the disk. The gas-to-dust ratio significantly lower than the canonical value (100 in 
the ISM) in the northern part indicates that the disk surface density (dust-fgas) is not high enough to collapse into a protoplanet by 
gravitational instability (Fukagawa et al. 2013). 

The radial distributions of dust and gas are fixed to Gaussian and power-law (with modifications at inner and outer radii) functions, 
respectively, in our modeling. The parameters that give a reasonable fit to observations are searched manually. Our model is one of 
the solutions but not unique, and we consider that the numbers listed in table 6 have at least a factor of ~ 2 uncertainties, especially 
for Eq. It is difficult to constrain the radial distributions more definately since the spatial resolution is limited. Still, we can argue 
that the gas distribution is more extended in the radial direction than dust and the gas-to-dust ratio is lower in the north than south. 

At a glance at figure 26, the gas distribution may be interpreted as “radially extended bump” and the dust particles are strongly 
accumulated in this bump region. Although the locations of the peak of the continuum emission and gas moment 0 profiles are 
different in the observed profiles (see figure 5), our model indicates that the peak of dust surface density profiles is within the region 
of gas surface density bump, which is between r ^ 100 AU and r ^ 250 AU. The peak of the gas emission at ~ 100 AU is partially 
due to the fact that the gas temperature is higher at inner radii. Interestingly, the bump in gas seems to be weaker in the south region 
than in the north. Such distribution of dust and gas seems to be, at least qualitatively, consistent with the picture of dust particles 
trapped in a pressure bump, or large-scale vortex, in gas (Pinilla et al. 2012; Birnstiel et al. 2013; Lyra and Lin 2013). We also note 
that recent observations by ATCA at 34 GHz indicate that large grains are concentrated in the northern part of the disk (Casassus 
et al. 2015). Vortices should be confined within a narrow radial range of the scale comparable with the disk thickness in this case. 
It is hinted that there might exist small-scale structures close to the peak of dust distribution, as discussed at the end of section 5.6. 
Future higher resolution observations may reveal the small-scale structures within the bright dust ring. Recently, Mittal and Chiang 
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Fig. 27. Optical depth profiles for continuum emission. Solid line is for the north profile and the dashed line is for the south profile. 


(2014) proposed that the global m — 1 mode might play an important role in producing radially extended lopsided structures. 
One interesting prediction of this model is that dust grains with different sizes accumulate at different locations within the disk. 
This may he addressed by higher resolution observations in multi-wavelengths and modeling effort including several dust species. 
Another interesting mechanism for the formation of dust-rich regions is the secular gravitational instability (Ward 2000; Youdin 
2011; Michikoshi et al. 2012; Takahashi and Inutsuka 2014). Takahashi and Inutsuka 2014 have done two-fluid analysis of gas and 
dust and shown that (multiple) ring-like structure with low gas-to-dust mass ratio can be created in self-gravitationally stable disks. 
Since initial gas mass should be relatively large in this case, significant gas dispersal should have occurred after the ring formation. 
This could be addressed, for example, by investigating the detailed kinematics of gas (e.g., Pontoppidan et al. (2011)). 


6.2 Optical Depth and the Total Amount of Dust and Gas 

6.2.1 Dust Optical Depth and Dust Properties 

We have seen that the dust particles are strongly concentrated in the north region. The radial profiles of the optical depth at the 
observed frequency (~ 340 GHz) along the line of sight is shown in figure 27, and it is clear that, in the northern part, the optical 
depth of dust emission reaches > 10. 

We have used the dust model with the maximum size of 1 mm. This is purely an assumption of this study, in order to estimate 
the minimum amount of dust grains needed to explain the bright thermal emission by using the dust model having the maximum (or 
at least large within the models that are considered to be reasonable) opacity at sub-mm wavelengths. In the case of dust model with 
Umax = 1 mm, the dust continuum emission is dominated by scattered light component if the disk is optically thick, which is the 
case for the northern part, because the dust scattering coefficient is much greater than the absorption coefficient (figure 6). Figure 
28 compares the radial profiles of dusf continuum emission between calculations with and without scattered light component. We 
have artificially sef the radiation energy density = 0 in equation (4) in the calculations to omit the scattered light component. In 
the northern part of the disk, the contribution from the scattered light component is a factor of ~ 5 larger than that from the thermal 
emission that has not experienced scattering. In the southern part, where the disk is optically thin, the contribution from the scattered 
light component is much smaller. 

There may be a variety of dust models that are able to reproduce observations. For example, the SED of HD 142527 can be 
reproduced by using irregularly-shaped micron-sized dust particles (Verhoeff et al. 2011). We have checked that the slope of SED 
in sub-mm range is consistent with observations in our model. Therefore, it is difficult to discriminate the dust size contributing 
most to sub-mm emission from current observations thus far. The scattering coefficient Ks of micron-sized dust particles at sub-mm 
range is expected to be much smaller than that used in our work. If the scattered light component can be observed exclusively, it is 
possible to discriminate the dust size. The polarization of dust continuum emission, for example, can be a good tracer of the dust 
size (Kataoka et al. 2015). 
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Fig. 28. Radial profiles of continuum emission calculated with (solid lines) and without (dashed lines) dust scattering. The left panel is for the north profile and 
the right panel is for the south profile. 

6.2.2 Gas Optical Depth 

We have used J — 3 —2 and C^®0 J = 3 — 2 emission to determine the gas structure. Figure 29 shows the radial profiles of 
the maximum optical depth of these lines (gas+dust optical depth, see also figure 30) for the reference disk models, ft is shown that 
the gas is optically thick at the peak of the lines. The maximum optical depth reaches an order of unity even for the south model of 
C^® 0 . 

Despite the fact that the lines are mostly optically thick, we have seen that the integrated intensity of gas emission becomes 
brighter as we increase the gas surface density (figures 11 and 12). This is because gas emission is nol entirely optically thick for 
all the velocity channels. When the amount of gas is increased, gas emissivity is increased and therefore, there are more velocity 
channels that contribute to the gas brightness at each given spatial location. Figure 30 shows the radial profiles of fhe optical depth 
of the line wings. We show the optical depth at the velocity channels that are different from those giving the maximum optical 
depth (line center) by 0.12 km/s, 0.24 km/s, and 0.36 km/s. For all of the calculations shown in figure 30, the optical depth of 
the continuum emission (which is derived using the first term of equation (10)) is subtracted from gas+dust optical depth (which 
is derived using all the terms using the right hand side of (10)) to show the gas-only optical depth^. The line emission becomes 
optically thin at 0.2 — 0.3 km/s away from the line center. Since the values of moment 2 shown in figure 23 and 24 are ~ 0.3 km/s 
at r ~ 150 AU, we consider that the line emission is optically thin when the velocity deviates from the line center by (only) half 
width of the line, which is close to the value of moment 2. Therefore, line wings can be used to trace the gas surface density even in 
the case that the emission at line center is optically thick. We therefore consider that we can reasonably constrain the amount of gas 
even when it is optically thick. However, we yet see that is more insensitive to the change of the amount of gas compared to 
which is more optically thin. Overall, we expect that the amount of gas has an uncertainty of a factor of ~ 2 — 3. For example, 
figure 11 indicates the models with Eq = 2.325 g cm“^ and Eo = 0.845 g cm“^ both explain the radial profiles of moment 0 in 
the north direction at r = 150 — 200 AU region reasonably well. To better constrain the amount of gas, we may need observations 
of lines with lower optical depth by, for example, using rarer isotopologues or using lower transition lines. The values of isotope 
ratio and the dust opacity are also important in accurately determining the amount of gas. Since there is fair amount of dust, the line 
emission is affected by the dust absorption, especially in the northern part of the disk. The scattering of gas lines by dust particles, 
which is not included in our model, can also be effective in determing the observed spatio-kinematic patterns. 

6.2.3 Total Amount of Gas and Dust 

We now discuss the total amount of dust and gas indicated from our modeling. We have two surface density distribution models, 
namely the north and south model, for each dust and gas. 

For dust mass, we obtain Mdust = 3.5 x 10“® Mq for the north profile and Mdust = 7 x 10“® Mq for the south profile when 
^ In other words, gas-only optical depth is calculated by taking into account only the second term of the right hand side of equation (10) 







26 


Publications of the Astronomical Society of Japan, (2014), Vol. 00, No. 0 


r[AU] 



0 100 200 300 400 

r[AU] 


Fig. 29. The line-of-sight optical depth of (solid line) and (dashed line). Top panel shows the optical depth of the northern part of the disk and the 
bottom panel shows that of the southern part. 
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Fig. 30. The line-of-sight optical depth of the north profile (left) and south profile (right). The upper panels show the optical depth of and the lower 
panels show that of C^®0. In each panel, the line-of-sight optical depth of the line center (the maximam optical depth of all the channels) is shown by black 
solid lines and that of the channels offset from the line center by 0.12 km/s, 0.24 km/s, 0.36 km/s is shown in red, green, and blue lines respectively. The 
optical depth of the continuum emission (black dashed lines) is subtracted for all of the gas optical depth. 
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integrated over the entire disk. Verhoeff et al. (2011) obtained the dust mass of 1 x 10~® Mq from their SED modeling. Considering 
that the bright arc of dust emission extends about ~ 1/3 of the whole azimuth, the value of dust mass we have obtained is similar to 
their results, despite the difference of grain properties. For gas mass, we obtain Mgaa = 1.8 x 10“^ Mq and Mgas = 5.7 x 10“® Mq 
for the reference north and south models, respectively, when integrated over the disk within 400 AU.® 

Although the total gas mass could have a factor of ~ 2 — 3 uncertainty as noted before, the overall gas-to-dust mass ratio may 
be < 10 — 30, which is smaller than the interstellar value of 100. The derived gas-to-dust mass ratio is likely to be an upper 
limit. Since the north region is optically thick, the dust emission is not very sensitive to the increase in the amount of dust (see 
appendix 2). We should also remember that the dust particles have relatively large opacity (and therefore emissivity) in sub-mm 
range (umax = 1 mm). If smaller grains are more abundant, we expect that the disk midplane, where dust emission originates, 
is colder while the disk surface, where optically thick gas emission originates, is warmer (Inoue et al. 2009). Consequently, even 
if small grains are more abundant, we expect brighter emission in gas while fainter in dust, resulting in smaller gas-to-dust ratio. 
We therefore propose that the dispersal of gas, or at least the depletion of CO gas, may have already taken place in the disk of 
HD 142527. However, it should be noted that detailed modeling with different dust grain properties is necessary to confirm this. 
Grain size and shape can strongly affect the observed emission properties (e.g., Min et al. 2003). 

Perez et al. (2015) derived the total gas mass of 0.1 Mq from their modeling effort. They used the model where the temperature 
is assumed to be 50 K at 50 AU with ® profile and the Doppler h parameter is assumed to be 50 m s“^. Their radial temperature 
profile is flat and the temperature is smaller than our model by a factor of ~ 1.7 at r = 100 AU. We have calculated how the 
assumptions of the temperature and the line width affect the results by setting those parameters to be the same as Perez et al. (2015) 
but the density is given by our reference model (equation (16) and table 6). The values of ^®CO J = 3 — 2 moment 0 is smaller 
by a factor of ~ 2 at r ~ 100 AU with this model. The impact of line width on moment 0 profile comes from the fact that there 
are less velocity channels that contribute to integrated intensity at each spatial position (see also the discussion in section 6.2.2). 
To compensate this by varying the amount of gas (parameter Eq in our model), we need to increase the value of Eq by a factor of 
~ 5 — 10, as seen in figure 11. Consequently, we consider that the difference of the total amount of gas comes from the different 
assumption of the temperature and the line width. The temperature profile used by Perez et al. (2015) comes from ^^CO J = 2 — 1 
observations, which may be smoothed by the relatively large beam. High resolution gas observations will play a decisive role in 
determining the distribution and total mass of gas component more accurately. 

6.3 Gas Bump Structures in Model 

The north radial profiles of the moment 0 of the reference model is bumpy at r ~ 170 — 200 AU, where dust emission is the brightest. 
As discussed in section 5.3, these bumps and troughs in gas emission arise because the continuum emission is very bright and comes 
from relatively narrow radial range. Since the line emission is partially hidden by dust, there is a significant effect on line emission 
when dust continuum is subtracted. 

There are at least two possible ways to obtain smoother moment 0 profiles. One possibility is to increase the amount of gas at 
the location of the peak of dust emission locally. In other words, there may be more detailed, small-scale (several tens of AU scale) 
structures in the radial direction than we have considered. As seen in figure 11, a factor of ~ 3 larger gas surface density locally at 
the peak of dust distribution may explain the observed values of moment 0 at this region. Since the dust emission is barely resolved 
in the radial direction, it may be difficult to see such small-scale variation of gas distribution. Higher angular resolution observations 
of gas is necessary to confirm this possibility. 

The other possibility is to consider more optically thin dust particles. Since the dust scattering opacity is ~ 10 times larger 
than the absorption opacity in our dust model, the total extinction coefficient of dust is dominated by the scattering opacity and 
is comparable with gas opacity (figure 30). Therefore, if the dust scattering coefficient is much smaller than considered in this 
paper, the gas emission can be brighter so the effects of continuum subtraction is more insignificant. Smaller dust particles may be 
one solution, but in this case, it is necessary to have larger amount of dust because the absorption coefficient (and therefore dust 
emissivity) is also small. 

A large (effective) emissivity and small scattering coefficient may be obtained simultaneously if we consider dust sedimentation. 
In our framework, the gas-to-dust ratio is constant in the vertical direction so large (~mm in size) dust particles reside even in the 
upper layer of the disk. However, such large particles may be sedimented to the disk midplane. The surface of the CO emission 
resides in the upper layer of the disk, and therefore there may only be small dust particles having small scattering coefficient around 

^ Inside 400 AU, the gas mass is dominated by the gas residing within r < r^ut although we have assumed that constant surface density in the radial direction 
for outer teneous gas. 
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the gas surface. The gas emission is not hidden by dust very much in this case and therefore the effects of the subtraction of the 
continuum emission is less significant. More sophisticated models that take into account the dust particle motion is necessary to 
verify this possibility. Also, the optical properties of the dust particles and its impact on observations should be carefully investigated. 
Detailed modeling in tandem with scattered light observations in near infrared (NIR; see also section 6.5), which is sensitive to the 
small dust particles in the upper layer of the disk, may be a key to verify this possibility. 

6.4 Discrepancy of Moment 1 at Inner Radii 

In the moment 1 profiles of the south gas models, we see that the reference model shows slower speed than observed at r < 100 AU 
for but faster at r < 50 AU for C^®0, as shown in figure 24. We have checked that this discrepancy is also present when 
the surface density profiles for the remaining gas inside the cavity is taken to be constant so the details of the radial profile of gas 
surface density do not alter the results. C^®0 data may be affected by the lower signal-to-noise ratio at high velocity channels, but 
the detection of at r ~ 50 — 100 AU region is robust. 

It is natural that the observed velocity should become smaller than Keplerian as one goes to inner radii, where the disk is barely 
spatially resolved in channel maps and the signal is weak. The gradient of line-of-sight velocity within the beam is large at inner 
radii and therefore emission at many velocity channels are averaged. Weak high velocity component may be discarded when the 
sensitivity is limited. These effects should properly be incorporated in our modeling since we convolve each channel map image by 
the Gaussian beam before calculating the model moment maps (see section 3.2). Therefore, the discrepancy between the observations 
and models might indicate that, in the southern region, the gas velocity is significantly faster than Kepler velocity at least for 
at inner radii. Rosenfeld et al. (2014) discussed that there may be a fast radial flow in the inner region of the disk of HD 142527 
based on and HCO*^ data taken by ALMA Cycle 0. Marino et al. (2015) suggested that the inner disk is significantly inclined 
relative to the outer disk based on the modeling of scattered light in the H-band. The discrepancy between the reference model and 
observations of for the southern region may also indicate the existence of such dramatic change of the inner disk structure. 
Observations with better spatial resolution are essential to clarify this point. 

6.5 Outer Tenuous Gas and Scattered Light in Near Infrared 

We have seen that there should be some amount of tenuous gas at r > 250 AU. The surface density of the tenuous gas is estimated 
to be ~ 7 X 10“® — 2 X 10“^ g/cm^. It is known that the HD 142527 disk shows an extended (~ 300 AU scale) scattered light 
emission in NIR observations (Fukagawa et al. 2006; Casassus et al. 2012; Canovas et al. 2013; Rodigas et al. 2014). Especially, the 
large-scale spiral structure appears from the south of the disk and extends towards the west. 

The extended scattered light emission in NIR may be connected to the existence of tenuous gas at outer radii. In NIR direct 
imaging observations, scattered light by small (typically, ~ 1 pm in size) dust grains at the disk surface is observed. Since such 
particles are well coupled with gas, we expect that small grains are distributed as gas is. The amount of tenuous gas (~ 0.02 g cm~^ 
at r ~ 300 AU in the south model) is probably enough to make the disk optically thick in NIR range, where dust opacity may be 
several tens of cm^ per unit gram of gas in the case of Umax ~ 1pm if gas-to-dust ratio is 100 in the outer r > 250 AU region (e.g., 
D’Alessio, Calvet and Hartmann 2001; Aikawa and Nomura 2006). Considering the detection limit of dust continuum, we expect 
that such (small) grains contained in the outer tenuous gas component are neither observable nor massive enough to contribute to 
the gas-to-dust ratio within 400 AU (see section 6.2.3). To further investigate the distribution of small grains, we need simultaneous 
modeling of scattered light and mm-emission, which is beyond the scope of this paper. 


7 Conclusion 

We have presented the models for gas and dust distribution for the disk around HD 142527 by using a series of axisymmetric models 
and by comparing the radial profiles of dust and gas emission. The schematic picture of the gas and dust distribution is shown in 
figure 31. Below, we list the main conclusions of our work. 

1. The dust particles are confined in a ring whose surface density profile is described by a Gaussian function in the radial directions 
with the width of Wd ~ 30 AU for both the north and south regions. The location of the ring is slightly different between the north 
{vd ^ 170 AU) and south profiles (r^ ~ 200 AU). To account for the strong azimuthal asymmetry of the continuum emission, 
the surface density of the dust ring is different by a factor of ~ 70 between the northern and southern part of the ring. Note that 
these results are based on the dust model with the maximum size of 1 mm, where the opacity at sub-mm wavelengths is large and 
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Fig. 31. Schematic picture of the disk around HD 142527. Green coiors indicate gas distribution and brown indicate dust distribution. The dust ring with the 
width of wj, 30 AU is iocated at r ~ 170 - 200 AU. The ring surface density is differ by a factor of ~ 70 in the northern and the southern part. The gas 
surface density is more radialiy extended. The gas distribution shows a bump-iike structure from ~ lOO AU to ~ 250 AU but there is remaining, tenuous gas 
component both inner and outer radii. 


therefore, we expect that the dust mass derived in our work is close to the lower limit. 

2. The gas distribution is azimuthally more symmetric than dust distribution. The amount of gas in the northern region is indicated 
to be ~ 3 times as large as that in the southern region. 

3. The gas radial surface density distribution is more radially extended than dust. The radial profiles of the moment maps of 

J = 3 — 2 and J = 3 — 2 emission can be described with the power-law profile between ~ 100 AU and ~ 250 AU, 
along with the remaining gas inside the cavity and with tenuous outer gas at r > 250 AU. Small-scale (~ 50 AU-scale) variation 
in the radial directions of gas distribution with a factor of ~ 3 is also indicated around the peak of dust distribution. 

4. Gas-to-dust ratio varies significantly within the disk. At the peak of the dust distribution, it is ~ 3 in the northern part of the disk 
while ~ 30 in the south. The gas-to-dust ratio over the whole disk is expected to be ~ 10 — 30 within 400 AU, indicating that gas 
depletion has already occurred in this system. 

We have used a series of axisymmetric disk models to do systematic parameter search while keeping the problem tractable. Since 
the CO isotopes used in this observation have turned out to be optically thick, observations of rarer isotopes or lower transitions 
are necessary to better constrain the amount of gas. More elaborate dust models will also be required to fully constrain the disk 
parameters and to account for observations at other wavelengths. 
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to improve the paper. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2011.0.00318.S. ALMA is a partner¬ 
ship of ESO (representing its member states), NSF (USA), and NINS (Japan), together with NRC (Canada) and NSC and ASIAA 
(Taiwan), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. 
This work is partially supported by JSPS KAKENHI Grant Numbers 23103004, 26800106 and 26400224. 


Appendix 1 Details on and results 

Figures 32, 33 and 34 show the channel maps of ^®CO from which the moment maps shown in figure 3 are created. Figures 35 and 
36 show the channel maps of C^®0 from which the moment maps shown in figure 4 are created. 
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As described in section 2.3, a constant (systemic) velocity of 3.7 km s“^ is found along PA = 71° ± 2°, which can be regarded 
as the direction of the minor axis of the system. The position-velocity (P-V) diagram along the major axis {PA — —19°; figure 37) 
was fitted by circular Keplerian motion for the emission detected above 5a (figure 37). The systemic velocity, position of the center 
of mass (the central star), and inclination relative to an observer are set as free parameters, whereas the central stellar mass is fixed 
in the range of 2.2 ± 0.3 Mq (Verhoeff et ah, 2011). The best-fitted parameters obtained by -minimization are in good agreement 
between and Using detected with a higher S/N, the system inclination angle is estimated as j = 26.5ti'Y degrees, 

where the uncertainty arises from the error in the stellar mass. Note that i is not large enough to yield reasonable constraints on both 
the stellar mass and the inclination (Simon et al. 2000). The systemic velocity is estimeted to be 3.72 ± 0.02 km s“^ in wlsr, and 
the obtained location of the center of mass matches that of the compact component of the continuum emission, which most likely 
represents the inner disk. Figure 37 also shows the curves for the Keplerian rotation with the parameters adopted in the modeling, 
i.e., i = 27°, M — 2.2 Mq and the systemic velocity of 3.7 km s“^ in ulsr- No significant deviation from the Keplerian motion 
was detected within the effective resolution of 0.2 km s“^ in our observations. 
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Fig. 32. Channel maps of (J = 3 - 2). The LSR velocity is shown in the top left corner of each panel, and the synthesized beam, 0"50 x o"42 with 
the major axis PA = 57.4°, is indicated by the ellipse in the bottom left corner of each panel. The white contours are the 5cr level, or 32 mjy beam^^. 
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Fig. 33. Channel maps of (J = 3 - 2), continued from figure 32. 
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Fig. 34. Channel maps of (J = 3 - 2), continued from figure 33. 
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Fig. 35. Channel maps of C^®0 (J = 3 - 2). The LSR velocity is shown in the top ieft corner of each panel, and the synthesized beam, 0"50 x 0''42 with 
the major axis PA = 60.6°, is indicated by the eilipse in the bottom left corner of each panel. The white contours are the So- level, or 41.5 mjy beam~^. 
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Fig. 36. Channel maps of (J = 3 - 2), continued from figure 35. 
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Fig. 37. The PV diagrams along the major axis (PA = 341°) in and (J = 3 - 2) lines. The contour spacing is 5a-, starting at the So- level. The 
dashed curves in each panel indicate the Keplerian rotation when M = 2.2Mq and the inclination angle i = 27°. 
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Appendix 2 Details of parameter search for dust distribution 

In this section, we present the results of parameter search for dust distribution models. We vary the parameters for dust distribution 
"rd and Wd and look for a set of parameters that best reproduces the dust continuum observations. Here, we exclusively show 
the results of the north region where dust emission is optically thick. For the south region, the disk is optically thin and the parameter 
search is straightforward. 

For each set of the parameters {T,d,o,rd,Wd), we calculate the dust continuum emission and the model images. The model images 
are convolved with the Gaussian beam with the size of the observations. The convolved radial profiles of the surface brightness is 
fitted by the Gaussian function given in equation (1). The derived parameters {Ip,ro,w) are compared with those derived from 
observations {Ip,obs,ro,obs,Wobs)- Tables 7, 8, and 9 show the results when ^d, and Wd are varied, respectively. 

The peak brightness of dust emission depends weakly on the assumed peak surface density. It is 1.17 and 1.25 ly/asec^ for 
Erf.o = 0.5 and 0.7 g/cm^, respectively, when other parameters are fixed at = 173 AU and Wd = 27 AU (the best-fit parameters; 
table 7). In other words, a factor of 1.4 increase in surface density results in only ~ 7 % in the surface brightness. The peak brightness 
depends weakly on and Wd as well. It is 1.24 and 1.19 ly/asec^ for Vd = 163 and 183 AU, respectively (table 8), meaning that 
12 % change in rd results in 4 % change in the peak brightness. In this case, ro also changes by ~ 20 AU so ro is more affected by 
the change in r^. When Wd is varied from 22 to 34 AU, the peak brightness changes from 1.12 to 1.28 Jy/asec^ (table 9), meaning 
that ~ 50 % change in the width of dust distribution results in ~ 14 % change in the peak brightness. In this case, the width of 
surface brightness changes by ~ 20 % so it is more affected than the peak brightness. 
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Model Parameters 

Model Results 

(Ed,o [g/cm^], ra [AU], Wd [AU]) 

{Ip [Jy/asec^], ro [AU], w [AU]) 

(0.5, 173, 27) 

(1.17, 152, 50) 

(0.6, 173, 27) 

(1.20, 153,50) 

(0.7, 173, 27) 

(1.25, 150, 50) 

Table 7. Results of model dust continuum emission of the north profile with different The 

{Jp,obs,ro,obs [AU],Waba [AC/]) = (1.2 Jy/asec^, 152 AU,51 AU). 


Model Parameters 

Model Results 

(Ed,o [g/cm^], Td [AU], Wd [AU]) 

{Ip [Jy/asec^], ro [AU], w [AU]) 

(0.6, 163, 27) 

(1.24, 141,51) 

(0.6, 173, 27) 

(1.20, 153,50) 

(0.6, 183, 27) 

(1.19, 161,51) 

Table 8. Results of model dust continuum emission of the north profile with different r^. The 

{Jp,obs,ro,obs [AU],'Wabs [AC/]) = (1.2 Jy/asec2,152 AU,51 AU). 


Model Parameters 

Model Results 

(Ed,o [g/cm^], Td [AU], Wd [AU]) 

{Ip [Jy/asec^], ro [AU], w [AU]) 

(0.6, 173, 22) 

(1.12, 154, 47) 

(0.6, 173, 27) 

(1.20, 153,50) 

(0.6, 173, 34) 

(1.28, 146, 56) 


The observed parameters are 


Table 9. Results of model dust continuum emission of the north profile with different ui^. The observed parameters are 

(Tp,o 6 s,r’ 0 ,o 6 s [AU],Wobs [AC/]) = (1.2 Jy/asec2,152 AU,51 AU). 
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